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Abstract 

We find pairs of solutions to a differential equation which is obtained as a special limit of 
a generalized spheroidal wave equation (this is also known as confluent Heun equation). 
One solution in each pair is given by a series of hypergeometric functions and converges 
for any finite value of the independent variable z, while the other is given by a series 
O ■ of modified Bessel functions and converges for \z\ > \zq\, where zq denotes a regular 

I singularity. For short, the preceding limit is called Ince's limit after luce who have used 

Oh' the same procedure to get the Mathieu equations from the Whittaker-Hill ones. We find 

^ as well that, when zq tends to zero, the Ince limit of the generalized spheroidal wave 
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equation turns out to be the Ince limit of a double-confluent Heun equation, for which 
solutions are provided. Finally, we show that the Schrodinger equation for inverse fourth 
and sixth-power potentials reduces to peculiar cases of the double-confluent Heun equation 
and its Ince's limit, respectively. 



o 

g: 1. Introduction 
in 

O ■ Firstly, we construct two linear differential equations whose solutions behave at infinity as the 

^ ! so-called subnormal Thome solutions, in contrast to solutions of a confluent and a double- 

^1 confluent Heun equations PP, from which the former equations are obtained by a limit process. 

I Secondly, we provide solutions which afford the expected asymptotic behavior for these equa- 
te ■ tions. Finally, we find that the Schrodinger equation with inverse fourth and sixth-power 
S ■ potentials reduces to particular instances of the double- confluent Heun equation and its Ince 
>- ! limit, respectively. 

I In the first place, let us introduce the two equations under consideration. Our starting point 

^ I is the generalized spheroidal wave equation (GSWE) in the form used by Leaver [2|, namely. 



cPU dU r 1 

z{z - zo)-^ + (5i + B2z)— + [^3 - 2r]u{z - Zo) + uo^z{z - zq)\ [/ = 0, (cu ^ 0) (1) 

where t] and uo are constants (notice that, if = and r] is fixed, the equation may 
be transformed into a hypergeometric equation). The points 2; = and z = zq are regular 
singularities with indices (0, 1 -|- Bi/zq) and (0, 1 — B2 — Bi/zq), respectively, while the infinity 
is an irregular singularity in which the behavior of U{z), inferred from the normal Thome 
solutions [S], is given by 

lim U(z) ~ e^'^^'z^'^'-^^-'/^l (2) 
2^00 

Since its parameters are not specified, the above GSWE is equivalent to the confluent Heun 
equation P], an equation that is more general than the original Wilson GSWE j^. Furthermore, 
as noted by Leaver, for 2:0 = we obtain a double-confluent Heun equation (DCHE) having 
five parameters, rather than four as in other contexts E] , namely, 

z"^ + (5i + B2Z)— + [53 - 2i^uz + u^z^] t/ = 0, (w 7^ 0, 5i ^ 0), (3) 
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where the singular points z = and z = oo are both irregular. For u = and/or Bi = (with 
T] fixed) this equation degenerates into confiuent hypergeometric equations (see Appendix A). 
At infinity, the behavior of U{z) is again given by 0, while at z = we find in the usual way 
that 

\imU(z) ~ 1 or \imU(z) ~ e^'^^z'^'^^ (4) 

z^O z^O 

The Leaver procedure also allows us to obtain solutions to the DCHE from solutions to the 
GSWE when zq goes to zero. The known Leaver-type solutions [2 E] are appropriate to solve, 
for instance, the Teukolsky equations for the extreme upper limit of the rotation parameter 
the time dependence of Dirac test fields in dust-dominated Friedmann- Robertson- Walker 
spacetimes and the Schrodinger equation with asymmetric double- Morse potentials [7j. They 
are suitable either for handling the Schrodinger equation with inverse fourth-power potentials, 
as we will see. 

Now, to get the equations we are interested in, the Levear limit is combined with a limit 
that Ince jH] had used to derive the Mathieu equation from the Whittaker-Hill equation. The 
Ince limit is obtained by taking 

a; ^ 0, r] ^ oo, such that 2riuj = —q, (5) 

where g is a constant. Thus, the Ince limit of the GSWE is 

d'^U dU 
z{z - ^o)^ + {Bi + B2z)— + [53 + q{z - z^)] U = 0, (q ^ 0). (6) 

This is a generalization of the Mathieu equation for, by setting 

zo = 1, Bi = -1/2, B2 = l, z = cos^ (au), W{u) = U{z), (7a) 
we obtain the equation 

d^W 

-— + [2q - 4^3 - 2q cos(2crM)] W = 0, (7b) 
du^ 

that is, the Mathieu equation if a = 1, and the modified Mathieu equation if a = i [H]. In 
fact, inserting zq = 1, Bi = —1/2 and B2 = 1 into Eq. (0), one recovers the algebraic 
Lindemann form for the Mathieu equation ^U] • Nevertheless, the trigonometric form ()7b|) with 
4i?3 = 2g — a is useful to verify that our solutions for the Ince limit of the GSWE give solutions 
already known for the Mathieu equation. 

On the other hand, the Ince limit of the DCHE - or Leaver limit of the Eq. (jH)) - is the 
equation 

(pTj du 
z^-^ + (5i + B2Z)— + (53 + qz)U = 0, (g 7^ 0, 5i ^ 0) (8) 

which degenerates into simpler equations if g = and/or i?i = (see Appendix A). Solutions 
are obtained for this special DCHE by taking the Leaver limit of solutions for the Eq. ©. By 
the way, we shall see that the Schrodinger equation for an inverse sixth-power potential is a 
particular case of Eq. (jSJ, as stated in the first paragraph. 

We emphasize that the Ince limits of the GSWE and DCHE, unlike the original GSWE and 
DCHE, require solutions behaving according to the subnormal Thome solutions jS], that is, 

lim U{z) ~ e±''v^z(i/^)'(^2/2)_ (g) 
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Despite this, our main mathematical issue consists in deriving pairs of series solutions to Eqs. 
(0) and (jSJ - having the behavior stipulated above at the singular points - from pairs of solutions 
to the GSWE. For this we shall again employ the Ince and Leaver limits. The solutions in each 
pair have the same series coefficients and these satisfy three-term recurrence relations. 

In section 2, a pair of solutions for the Ince limit of the GSWE is obtained by taking the 
Ince limit ((51) of a known pair of solutions for the GSWE. One solution is given by an expansion 
in series of hypegeometric functions and converges for any finite z; the other solution is given 
by an expansion in series of modified Bessel functions and converges for \z\ > \zq\. Other pairs 
are generated by using transformation rules. These rules result from variable substitutions that 
preserve the form of the differential equations but modify their parameters and/or arguments. 

In section 3, we find pairs of solutions for the Ince limit of the DCHE by taking the Leaver 
limit (2:0 0) of solutions for the Ince limit of the GSWE. Solutions in series of irregular 
confiuent hypergeometric functions result from expansions in series of hypergeometric functions 
and converge for any finite z. The other solution in each pair is given by a series of modified 
Bessel functions and converges for \z\ > 0. 

In both of these sections we deal with solutions with and without a phase parameter u. In 
general, this u is introduced in order to assure the convergence of the series when there is no 
free constant in the differential equation, as in some scattering problems or in equations where 
^ is a variable related to the time [2,7. Solutions with a phase parameter are two-sided in 
the sense that the summation index n runs from — cxo to 00. However, if there is an arbitrary 
parameter in the equation, we can truncate the series by requiring that n > 0. In this manner, 
we obtain u in terms of parameters of the differential equation. 

In section 4, we show that the Schrodinger equation with inverse fourth and sixth-power 
potentials in fact leads to the DCHE and its Ince limit. Some additional considerations are 
provided in section 5, while in Appendix A we discuss the degenerate cases of the DCHEs, in 
Appendix B we present an alternative derivation of the expansions in Bessel functions, and in 
Appendix C we rewrite the Leaver-type solutions for the DCHE in a form appropriate to solve 
the Schrodinger equation with an inverse fourth-power potential. 

2. Ince's limits for the generalized spheroidal wave equa- 
tion 

In this section we use transformation rules that permit us to generate new solutions from a 
given solution for the Ince limit of the GSWE. The rules Ti, T2 and T3 below are derived 
from the ones valid for the GSWE [7j and can be checked by substitution of variables. If 
U{z) = U{Bi, B2, B3] zq, q; z) denotes one solution for Eq. (jH)), the effects of these rules are as 
follows 

m{z) = z'+^^/^<^U{Cu C2, C3; zo, q; z), zo^ 0, 

T^Uiz) = {z- ^o)'-^^-^^/^"f/(fii, D2, D3; Zo, q; z), zo + 0, (10a) 
T3f/(z) = U{-Bx - B2Zo,B2,B3 - qzo;zo, -q; Zq - z), 

where 

C, = -B,-2zo, ^2 = 2 + ^2 + ^, C3 = i?3+(H-f)(52 + f), 

\ ^oJ\ zoJ ^^^^^ 

D2 = 2-B2-'-^, ^3 = i?3 + f (f + 52-1). 

We use only Ti and T2. The rule T3 exchange the position of the regular singular points 
z = Zo ^ z = and may be used to get an alternative representation for the solutions, but 
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these are not proper for getting the hmit Zq ^ 0. 

In section 2.1 we derive two pairs of solutions for the Ince hmit of the GSWE - denoted 
by (t/j°,f/°f), i = 1,2 - with a phase parameter u. The superscript 'zero' indicates that the 
series converges in any finite part of the complex plane, while the superscript 'infinity' indicates 
convergence for \z\ > \zo\. The second pair of solutions results from the first by means of 
the rule T2. In section 2.2, we truncate these series by taking n > and obtain four pairs of 
solutions without phase parameter. 



2.1. Solutions with a phase parameter 

Denoting by 6^ the series coefficients of a solution, their recurrence relations will have the 
general form 

anK+l + PnK + InK-l = 0, (-CX) < Tl < Oo) (11a) 

where a„, 7„ and 6„ depend on a phase parameter u which may be determined from a 
characteristic equation given as a sum of two infinite continued fractions, namely, 

_ a_i7o q;-27-i «-37-2 , ao7i «i72 «273 

For a specific pair of solutions we add a superscript in each of these quantities. 

The first pair of solutions for the Ince limit of the GSWE comes from the following pair of 
solutions of the GSWE |7j 

n=-oo \ ^ Z Zo Zo 

(12a) 

U^{z) = e*^^2i-(^2/2) J2 bll\-2iujz)'^^''^{n + v + I + ir],2n + 2v + 2] -2iujz), 

n=— oo 

where F{a,b;c;y) and '^{a,h]y) denote, respectively, the hypergeometric functions and the 
irregular confiuent hypergeometric functions I12j. The solution Ul^, converges for any finite 
z, whereas converges for | z |>| |. In the recurrence relations ()lla|) for h\^^ we have 



2{n+v+l)in+v+ 



riLOZQ\ 



(12b) 



= -B, - ^uz, - (n + + 1 - f )(n + + f ) 

^ = -^^^0 2(nVl)(n+.) " 

Note that u cannot be integer or half-integer in order to avoid vanishing denominators in the 
coefficients of the recurrence relations. Moreover, for an integer or half-integer u, we would 
have two equal hypergeometric or confiuent hypergeometric functions (for different values of 
n), contrary to the hypothesis that all the terms of the series are independent. On the other 
hand, the hypergeometric functions are not defined if B2 + (Bi/zq) is zero or a negative integer. 
Nonetheless, a transformation rule supplies another solution which is valid for these values of 
B2 + (Bi/zo). 

The three-term recurrence relations ()lla|l constitute a infinite system of homogeneous linear 
equations for which nontrivial solutions for the coefficients 6„ demand that the determinant 
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of respective tridiagonal matrix vanishes. Equivalently, the characteristic equation must be 
satisfied and this is a condition necessary also to assure the convergence of the series by means of 
a Poincare-Perron theorem JH]- However, there are two possibihties to satisfy this requirement. 

On the one hand, if there is some free constant in the differential equation, that constant 
must be determined so that the characteristic equation is fulfilled for the admissible values of v 
(that is, neither integer nor half-integer). In this case, the freedom of choosing v may be used in 
two different ways: (i) to obtain two-sided solutions (—00 < n < 00) by ascribing appropriate 
values for i/, or (ii) to obtain one-sided solutions by choosing p such that n > 0. At the end 
of the present section, we use the first alternative to rederive some Poole's solutions 113] 
for the Mathieu equation, having period 27rm, where m is any integer equal or greater than 2. 
In section 2.2, we use the second alternative for the general case. These latter solutions afford 
solutions with period vr or 27r for the Mathieu equation, in contrast to the solutions obtained 
in the first alternative. 

On the other hand, if there is no arbitrary parameter in the differential equation, the 
parameter v takes the role of free parameter in the sense that it must be adjusted to ensure the 
validity of the characteristic equation and, consequently, the convergence of the series. By this 
reason, v is also called characteristic index or parameter ^H] ■ Examples of equations requiring 
a phase parameter are discussed in section 4. 

From the above pair of solutions for the GSWE, by using the Ince limit (0), we readily find 
the solution Ui^{z) written in the first pair below. To get the Ince limit of the solution U^{z)^ 
we define Cn as 

6« = {ir^Y+''T{iri - n - i/)c„. 

This imply that 

n=-oo \ 

where q = —2riuj. The recurrence relations for c„ are 

OinCn+l + (3^n^Cn + InCn^l = 0, (-00 < n < Oo) 

with 

ir] _ ir]-u-l 

trj — n — u ~ 1 17] 



On the other hand, we have 

lim [r(a + 1 - c)^(a, c; x/a)] = 2x^^"^l'^K^A2^) (13) 

a— +00 

where K\{^ denotes the modified Bessel function of the second kind ^12j whose definition in 
terms of irregular confiuent hypergeometric functions is il2j 

Ky^iX) = K^xiO = 6-^20'^ (a + i, 2A + 1; 2^) . (14) 
Then, using ()13|) we find that for ir] oo {n fixed and q =constant) 



r(^r/ -n- v){-qzY^-m [n + u + 1 + t7],2n + 2u + 2; -|) ^ 2{-qzy/'K2n+2u+ii±2t^), 
lima„ — >■ lima^^-*, lim7„ lim7,^^) =^ limc„ — > lim6^^^ 
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Using these results, we find the Ince hmit of U^, written in the first pair below. Although this 
is a formal derivation, the solution may be checked by inserting it into Eq. © (see Appendix 
B). In addition, from the relation 

we see that the expansions in series of Bessel functions have the behavior given by 

jimf/-(^)~e^^^v^z(V^)-(^^/2), < arg(±2.y^) < (j = l,2) 

in accordance with Eq. Q . The second pair of solutions follows from the first one through the 
rule T2, as mentioned before. Moreover, solutions for the Mathieu equation are obtained by 
using Eqs. ()7aj) and by noting that in this case the hypergeometric functions can be rewritten 
in terms of trigonometric functions. 



First pair 



UU^)= E -n-.-l,n + .+ f;5, + ^;l-f 

n=-oo ^ ^ Zo Zq, 

00 



(16a) 



where in the recurrence relations (lllal 



(„+,+2-^)(n+.+l-^-^ 
(n+y+l)(n+i/+|) 



= qZo- 



(n+!^-i)(n+;^) 



If B2 + {Bi/zq) is zero or a negative integer we have the solution U21, instead of 
For the Mathieu equation we use Eqs. (fTaj) and the formula [TT] 

F —a, a; (1/2); sin^ (ctm) = cos(2ao"M). 

Thence, we obtain even solutions with respect to u, namely, 

00 

W^^{u) = J2 cos[(2ra + 2u + l)au\, \ cos(a-M)| < 00, 

n=—oo 
00 

WS^iu) = E b^n^K2n+2u+i [±2i^ cos{au)] , I cos((tm)| > 1, 

?1=— 00 

with the simplified recurrence relations {a = 2q — AB^) 

qh^^l^ + \{2n + 2v + if - a] 6« + ql^^l^ = 0. 



(16b) 



(17a) 



(17b) 
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Second pair 

B S ^ 



X 



F [-n-u- f + l,n + iy + 2- f;2- B2 - 



:i8a) 



^ ^0 z^o^ 2 2 



where 



a(2) = g2:o 



7i^^ = qzo 



(n+i/+l)(n+i/+|) 

(n+.+l^^)(„+.-^-£L 
(n+i/-i)(n+;/) 



;i8b) 



in the recurrence relations for 6(f). If B2 + (^i/^o) 

is a positive integer equal or greater 
than 2 we have the solution U^j^ instead of 1/2^. Note that, in writing the solution 1/2^, we have 
used the relation 



F(a, b; c; y) = (1 - yf-'^-'Fic - a, c - b; c;y). 

For the Mathieu equation we use the relation JT] 

, - . 2/ N \ sin[(2a - l)au] 

F a, 1 - a; -; sin^fau) ^ - 



(19) 



(2a — 1) sin(o"M) 

and, in addition, define c„ as b^^ = (2n + 2t/ + l)c„. So, we find that the recurrence relations 
for c„ become identical to the ones for b^\ giving the odd solutions 



Wl{u)= E &i'^sin[(2n + 2z/ + l)aM], 

n=— 00 

00 

W^{u) = tan (au) ^ (2« + 2^/ + 1) bll^K2n+2u+i [±2i^cos{au)] 



(20) 



where | cos(o"n)| < 00 and | cos(crM)| > 1, respectively. 

As we have explained earlier, if there is a free parameter in the differential equation, it 
is possible to satisfy the characteristic equation for any noninteger or half-integer u. We use 
this fact to rederive some Poole's solutions ITH] to the Mathieu equation. For this, in the 
previous iy{[,(M) and iy{|,(M) we take 



2z/ + 1 = //m, a 



(21) 



where / and m are integers prime to one another, / < m. Then, we find the two-sided Poole 
solutions Wfi^u) and Wf{u) given by 



W[{u) = E cos 



2n + — ]u 
m 



sm 



2n + — ]u 
m / 



(22) 



The first is even with respect to u and the second is odd, and both of them have period 27rm, 
m > 1. Since they have the same series coefficients, we can combine them to find another Poole 
solution, that is. 



W''{u)= E ^i'^exp 
?i=— 00 
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i \2n -\ ] u 

\ ml 



(23) 



Furthermore, for an arbitrary v we find 



oo 

W[u) = (2^ + 2z/ + 1) m] , (24) 

n=— oo 

which is also a solution already known in the literature ITK] . 



2.2. Solutions without phase parameter 

Now we truncate the solutions obtained in section 2.1 by taking n > 0. This gives v in terms 
of some parameters of the differential equation. The resulting solutions are convergent only if 
there is a free parameter to be determined from the characteristic equation. 

This truncation reverses the procedure by which the solution Uij^{z) for the GSWE, given 
in Eq. ()12a|l . was obtained. Indeed, that solution was constructed jTHj as a generalization of 
an one-sided series of Jacobi polynomials, constructed by Fackerell and Grossman (TH] to solve 
the angular Teukolsky equations of the relativistic astrophysics. Despite this, the truncated 
solutions found in 7J are more general than the Fackerell-Crossman ones because no particular 
values are attached to the parameters of the GSWE and also because the truncation was 
extended to the Leaver expansion U^{z). I addition, these one-sided series are suitable to get 
solutions in finite series, the so called quasi-polynomial solutions. In effect, a solution whose 
coefficients 6„ obey recurrence relations as 

OinK+l + l3nK + InK-l =0, 72 > 0, 6_i = 



becomes a quasi-polynomial solution with < n < — 1 whenever 7Ar = for some n = N 

m ■ 



For the truncated solutions - denoted by (f/°, f/°°), i = 1, 2, 3, 4 - the recurrence relations 



and the characteristic equations have one of the three forms written below. The first case 
(a;_i = 0) is the general one and the others (a-i 7^ 0) may occur only for special cases. 



ao&i + Pobo = 0, 

a„6„+i + PnK + Inhn-l = {u > 1), 



/5o 



«07l <^ll2 «273 

A- (32- /3s- 



(25) 



aobi + Pobo = 0, 

ai&2 + Pibi + + 7i] bo = 0, 

a„6„+i + Pnbn + Inhn-l = (n > 2) 



/5o 



ao + 7i] ai72 0273 



/32- /53- 



(26) 



aobi + [/5o + 60 = 0, 

anbn+l + (3nbn + 7n&n-l = (w > 1) 



+ a. 



ao7i Qi72 0^273 
A- Ps- 



(27) 



Note that we have n > — 1 in a„, n > in /5„ and > 1 in 7„. 

These forms for the recurrence relations are the same that appear in truncation of the 
expansions (jl2aj) for the GSWE j7]. As a matter of fact, the solutions of the present section 
are the Ince limit of solutions for the GSWE given in section 3 of Ref. [7j. However, in order to 
illustrate how these recurrence relations are obtained, we insert the solution given in ()16a|) 
into the Ince limit of the GSWE. Then, for n > 0, from Eq. ()B4|) we find that 



oo 

E 

n=0 



2n+2i^-l 



(0 



n=0 



n=0 



2n+2u+3 



(0 = 0. 
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Setting m = n — 1, m = n and m = n + 1 in the first, second and third terms, respectively, this 
equation becomes 



a_i6o^2.-i(0 + K&i + Pobo] i^2.+i(0 + [c^ih + Pih + 7i&o] ^2.+3(0 + 

oo 

[anbn+1 + Pnbn + 7n&n-l] K2n+2u+l{0 = 0) (2^ 



n=2 



where we have dropped the upper suffixes. Therefore, if we can choose u so that q;_i = 0, we 
find the ffist set of recurrence relations. However, notice that 

g.o(. + l-f)(.-^-f) ^ . r . = f -lforB,^l,2; 

2.(.+ l/2) ' l- = t + f forf + f ^0,1- 

Hence we see that there are two possible choices for u and, for each of them we have two cases 
in which a_i may differ from zero. Let us consider only the case u = {B2/2) — 1. Then, for 
the exceptional case B2 = 1 {u = —1/2), we find K2U-1 = K2U+3 = K2 (since Kx = K_\) and 
therefore the Bessel functions in the ffist and third terms of Eq. (j^Hj) are equal, giving the 
recurrence relations (j26|) . Similarly, ii B2 = 2 {u = 0) we find K2u^i = K2U+1 = Ki and this 
leads to the recurrence relations (|27|). In this manner we obtain the ffist pair given below. The 
remaining can be derived from this by using the transformations rules Ti and T2 as 

The condition on each pair is imposed in order to assure that the special functions are inde- 
pendent in both solutions; it guarantees either that there is no vanishing denominator in the 
recurrence relations. Furthermore, we have additional restrictions on the parameters of the 
solutions Uf . Thus, if B2 + [Bi/zq) is zero or a negative integer, the hypergeometric functions 
are not defined in and f/g t>ut are defined in f/3 and f/4 , and vice- versa. The results for the 
Mathieu equations are already known [^j, but note that the recurrence relations for this case 
come from the three Eqs. (j25H27|) above. 



First pair: B2 7^ 0, —1, —2, ■ ■ -. This ffist pair corresponds to u = {B2/2) — 1 in (t/i^,, U[ 

^?W = E^i'^^ + ^2- i;i?2 + —;!-- 

n=0 ^ ^0 Zq 

00 

Ur{z) = Y: b^^^K2n+B2-l (±2^v^) , 



(29a) 



n=0 



with the following coefficients 



/3<') = 4B, - 2,.„ + 4n („ + - 1) - , (2*) 

^" = („+i^_|)(n+i^_i) 

in the recurrence relations for the b^^\ namely: Eqs. (j2SD if ^2 ^ 1, 2; Eqs. ^ ii B2 = 1; Eqs. 
(EZD if = 2. 
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For the Mathieu equation we find solutions 

oo 

W^iu) = X! cos{2nau), \ cos(cr-u)| < oo, 

n=0 



oo 



(30a) 



W^{u) = J2 b^nK2n [±2i^cos{au)] , I cos(aM)| > 1, 



n=0 



with the simphfied recurrence relations 

qb^^^ - ab^o'> = 0, qb^2^ + [4 - a] b^^'> + 2qb\^^ = 0, 

qb^^l + [An' - a] 6« + qb^^i = (n > 2). 
These solutions are even with respect to u and, for o" = 1, the solution W^{u) has period tt 



(30b) 



Second pair: {B2/2) + (Bi/zq) 7^ —1, —3/2, —2, —5/2 ■ ■ -. This pair of solutions can also be 
obtained by taking u = {B2/2) + [Bi/zq) in {U^^, U^^). 



Ul^z) = ^i+(^iAo) ^ 5(2)^ f_n^ ri + l + B2 + ^;B2 + ^;l-- 

\ Zq Zq Zq 



n=0 

00 



(31a) 



U^{Z) = Y: &i'^i^2n+l+B.+(2B,M,) (±2^v^) 

n=0 



where 



(2) 



qzo{ 



/3f = 4^3 - 20^0 + 4:(n + l + ^) (n + B2 + ^) - , , 



' \'^^oJ (31b) 



(2) 
7n 
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in the recurrence relations for b'^^: Eqs. if (^2/2) + (Bi/zq) ^ 0,-1/2; Eqs. ^ if 

(S2/2) + (Si/^o) = -1/2; Eqs. (EID if (^2/2) + {B^jz^) = 0. 
For the Mathieu equation we again have even solutions 



COS^cr-UJl < 00, 



l^2°(w) = E^n^ cos[(2n + l)H, 

n=0 

00 

W^iu) = E b^^K2n+i [±2t^cosiau)] , | cos(an)| > 1, 



(32a) 



n=0 

with the recurrence relations 

^(2) , U , 1 ^l/,(2) 



qbY' + [q + l-a] b\;> = 0, 



qb^^l + {2n + lf-a 6f + ^^^(f)^ = q (n > 1). 



(2) 



(32b) 



If (T = 1 the solution (m) has period 27r. 
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(33a) 



Third pair: B2 7^ 4, 5, 6, ■ ■ -. This corresponds to u = 1 - {B2/2) in {U^^, U^). 

B B / R r 

Ul{z) = {z- z.f-'^'-^z'^^ ( -n, n + 3 - 2 - 52 - — ; 1 - - 

n=0 ^ ^0 Zq 

n=0 

with the coefficients 

= 4i?3 - 2gzo + 4(n + l)(n + 2 - B^) - ^J,l^)(ll2-!^) ' ^^^^^ 

in the recurrence relations: Eqs. ^ if B2 ^ 2, 3; Eqs. ^ if B2 = 3; Eqs. (|17|) if ^2 = 2. 

For the Mathieu equation we redefine the coefficients b!f^^ as b!f^^ — * (2n + 2)6^^\ Then we 
find the odd solutions 



Wf^{u) = sin[(2n + 2)au], \ cos(aM)| < cx), 

n=0 

34a 

00 ^ ' 

2„+2 [±2ii/g cos(cru)] , I cos(crM)| > 1, 



n=0 

with the recurrence relations 

qbf^ + [4 - a] b'i^ = 0, 

(34b) 

gfeSi + [4 (n + 1)' - a] 6(f) + g^i^^li = 0, (n > 1). 
For a = 1 the solution (m) has period vr. 

Fourth pair: {B2/2) + (i?i/2;o) 7^ 1, 3/2, 2, 5/2 ■ ■ -. This can also be obtained by setting 
u = -{B2/2)-{B,/zo) in {U^^U^,) 



U! = {z-zot'''-^f:b^:^F(-n,n+l-B2-'^;2-B2-^;l-- 

\ Zq Zq Zq 



n=0 



(35a) 



n=0 

with coefficients 



qzQ (n+l)^n--^^ 



/jW = 4i?3 - 2gzo + 4 (n - f ) - i?2 + 1 - f ) - ^^ .^^^^^ /^.U - (35b) 

(4) ^ !!lH:f!li)H!llL 



7. 



"2 2 M " 2 ^0 
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in the recurrence relations: Eqs. 1^ if (52/2) + (Bi/zq) ^ 0,1/2; Eqs. 1^^ if (^2/2) + 
(Bi/zo) = 1/2; Eqs. (EH) if (^2/2) + (Bi/zo) = 0. 

For the Mathieu equation we redefine 6„(4) according to 6^^-* — * {2n + 1)6^^-* and find the 
odd solutions 



W^{u) = J2 ^ sin[(2n + 1) 



cos aw < 00, 



n=0 



(36a) 



W^{u) = tan (au) ^ {2n + 1) foi^^Kan+i [±2zv^cos(aM)] , | cos(crM)| > 1, 

n=0 



with the recurrence relations 

r(4) 



(4) 



564 + [1 — g — a] 6| 
rfli + \{2n + 1)^ - a] foW + g6l^2i = (n > 1). 



(36b) 



Now, for cr = 1, W^{u) has period 27r. 



3. Ince's limits for the double-confluent Heun equation 

As in the case of the Ince limit of the GSWE, we have found no solution in the literature for the 
Ince limit of the DCHE. The solutions below are obtained by taking the limit zq (Leaver 
limit) of the solutions given in section 2 for the Ince limit of the GSWE. For this we use the 
formulas |T2j 

jinn F ^a, 6; c; 1 - ^ j = ?/"^(a, a+l-b;y), (37a) 



lim 1 + ^ = ^ lim 1 - - = e^^/^ (37b) 

V a J zo-^o V z J 

Actually, it is not necessary to use the second equation above, since we can get one pair of 
solutions as the limit of the first pair of section 2.1 and, then, generate the other pair by means 
of the transformation rule 

rU{z) = e'''/'z^-''^U{-Bi,4: - ^2, ^3 + 2 - B2] q; z), (38) 

where U{z) = U{Bi, B2, B-y, q; z) denotes known solutions of Eq. (jH)). On the other hand, to 
check that the solutions U^{z) exhibit the behavior given in Eq. (0)) when z — > 0, we may use 
the relation ^2] 

\im ^{a,b;y)r.y-^[l + Oi\y\-\ - < argy < — . (39) 
\y\-*oo Z Z 

3.1. Solutions with a phase parameter 

For the solution Ui^, of section 2.1, we find that the limit of the hypergeometric functions when 
Zq tends to zero, B2 and Bi being fixed {c = B2 + Bi/ zq ^ 00), is given by 

^ / B2 B2 ^ Bi z 
lim F [n + u -\ , —n — u — 1 -\ ; B2 H ; 1 

20^0 V 2 2 Zq Zq 

oc z-^-^ ^ (n + u + ^,2n + 2u + 2;^] . 
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Then, considering also the solution and the limits for the coefficients in the recurrence 
relations, we get the ffist pair of solutions with a phase parameter u (different of integer or 
half-integer). The rule r leads to the second pair. 



First pair 



n=—co 
oo 



z J \ 2 z 



(40a) 



f/-(z) = Y: &i^)i^2n+2.+l (±2zy^) , 

n=— oo 

where in the recurrence relations piaj] 

^ gBi(n+.+2-^) 

/?(^) = 4i?3 + 4 (n + . + 1 - f ) (n + . + f ) - > (40b) 

(n+z/-i)(n+i/) ■ 

Second pair 

UU^) = e^^/^z-'^-^ £ 6(2)f-:?i)%L + z/ + 2-|^,2n + 2z/ + 2;-:^ 



2 1—7 



oo 



(41a) 



n=— oo 

where 

«r = Pn = In = , (41b) 

in the recurrence relations ()lla|l for 6^^^^ 

3.2. Solutions without phase parameter 

These solutions may be derived by truncating the solutions of section 3.1. In this case, we see 
that there is only one choice for v in each pair. Alternatively, the solutions can be found by 
applying the Leaver procedure to the ffist and third pairs of section 2.2. 

First pair: B2 ^ 0, -1, -2, ■ ■ -. This corresponds to v = {B2/2) - 1 in (f/°^, U^). 



°° / R \ " / R 
U^{z) = z'~^'^ Y^b^nH — ) ^ ( ^ + ^2 - 1, 2n + B2; ^ 



n=Q 



Z / \ Z 



(42a) 



n=0 

with the following coefficients 

K?)R?Ti)' 



= AB, + 4n (n + 5, - 1) - ^^.t^^.^y (42b) 



(1) _ gBi(n+B2-2) 
7n - („+^_3)(„+^_i) 
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in the recurrence relations for the b^^^: Eqs. ^ if B2 ^ 1, 2; Eqs. ^ if B2 = 1; Eqs. ^ if 
B2 = 2. 

Second pair: B2 7^ 4, 5, 6, ■ ■ -. It corresponds to u = 1 - {B2/2) in {U^^, U^) but can also be 
obtained from the first pair via the rule r. 

U^[z) = e^i/^^-i bi? l-^) ^ U + 3 - 52, 2r2 + 4 - B2; 



n=0 



n=0 



where 



(3) ^ gBi (n+1) 



(43a) 



V-+2-^)(n+|-^)' 

/3(3) = 4S3 + 4(n + l)(n + 2 - B2) - (^.^ff^l;!^) ^ (43b) 

^(3) _ gBl («+2-B2) 

7n (n+i-ia)(„+i„ia)- 

in the recurrence relations for h^^^. Eqs. (|23) if S2 ^ 2,3; Eqs. (1211) if ^2 = 3; Eqs. (j27j) if 
^2 = 2. 



4. Potential applications 

As we have mentioned, the Schrodinger equation with inverse fourth and sixth-power potentials 
can be reduced, respectively, to the double-confluent Heun equation Q and its Ince limit (jH)). 
Singular potentials like these have appeared in the description of intermolecular forces [SB] and 
in the scattering of ions by polarizable atoms. For the sake of illustration, we consider the last 
problem. 

Before discussing these examples, let us present the so called normal forms of the DCHE, 
that is, the forms in which there is no first-order derivative terms in the differential equations. 
The general procedure for this, consists in writing the equation as 

+ p{z)— + q{z)U = 0. 



dz"^ dz 



Then, the substitution 



U{z) = F{z) exp ^— - j p{z)dz 



gives a first normal form, namely. 

From this, other normal forms are obtained by the transformations 

fdh 

z = h{d), F{z) = G{^) 



d§ 



which yield 



d'^G (dhV Id^h ,dh 3 fd'^h ,dh^ ^ 



„.^mG=o, m^nm]i- +__/___ _/_ 
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By employing this procedure, Lemieux and Bose j2I] have derived several normal forms for 
the general Heun equation and its confluent cases, excepting the triconfluent equation. These 
forms are useful to recognize whether a given equation belongs to the Heun class. Nevertheless, 
to find the solutions for the equation, we have to come back to the form for which the solutions 
were established, as below. The three Lemieux-Bose normal forms for the DCHE, together with 
the transformations of variables, are the following: 



U{z) 



~B2/2^Bi/{2z) p/ 



dz 



a p I 

J, 2 "T 



0; 



(44) 



dp 



J„2 "T 



G = 0; 



(45) 



e^", f/(2) =//(«) exp[^A(l-52)M + ^e-H ^ //(m) = z(^2-i)/2g-i?i/{2^)^(^^) 



d^H 



+ A2 



5, 



I--B2 r _ :^p-2A« _ u (B2 
2 I A ^ -"112 



1 ) e"^" - 2r/u;e^" + u'-'e 



Xu 



,2„2Au 



if = 0, 



(46) 



where A is a constant at our disposal, for example, A = 1 or A = Note that, since these 
transformations involve neither 77 nor u, their Ince limits are obtained by putting uj"^ = and 
2riuj = —q. 

Now we proceed with the scattering problem. The radial part R{r) = xi'f')/^ of the wave 
function for the Schrodinger equation in three dimensions, for a particle with mass /i and energy 
E, is 



(fx{r) 



+ 



r2 h 



X{r) = 0, 



(47) 



where fc^ = 2pE/h^ , I is the angular momentum and V{r) is the potential. Now, according to 
Kleinman, Hahn and Spruch 22 , for the interaction of a light particle of charge e' with a fixed 
atom of charge Ze containing z' electrons, we have 



V{r) 



{Z - z')ee' 



2r4 



- (^2 - QaoP'i) 



-,'\2 



2^6 



(4^ 



where r is the distance from the incident ion to the atom, = /i^/(/ie^) is the Bohr radius, 
a[ and a'2 are, respectively, the electric dipole and quadrupole polarizabilities of the atom and 
I3[ is a parameter resulting from a nonadiabatic correction {a'l, 0L2 and are constants which 
describe the properties of the target only). For this potential, the Schrodinger equation becomes 



+ 



2fi{Z 



z'W 



+ 



pa\{e 



l\2 



p {a'2 — 6ao/?i) (e 



l\2 



x = o. 



(49) 



Therefore, for neutral targets {Z = z') this is a particular case of the Ince limit of the DCHE, 



as we see from Eq. (jl^j) with 



0, 2riuj = —q and z = p^ = r"^. On the hand, if the inverse 



sixth-power term vanishes (og = 6ao/5i), this radial Schrodinger equation is a particular case 
of the DCHE as seen from Eq. with B2 = 2, for neutral or ionized targets. In both cases 
the energy of the incident particle (/c^) is given and, consequently, there is no free parameter 
in these equations since the other constants are also fixed. Then, convergent solutions require 
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a phase parameter i/, analogously to the scattering by the field of an electric dipole 0. To 
obtain the radial dependence R{r) we must convert Eq. ()49|) into the DCHE Q and its limit 
(jH)). Below we discuss only the asymptotic behaviors of the solutions for the each case. By this 
reason, we do not write the recurrence relations for the coefficients. 

Potential with inverse fourth and sixth-power terms. Eq. ()45|1 suggests the substitutions 
which transform the Schrodinger equation (jl^ into 



,(fU 



dU 



z--— + {Bi + B2z)— + 
az'^ az 



B2 

2 



^2 _ 3 
2 4 



'-^+'-z-^iZ-z')V-z 
4 4 2^ 



U = 0. 



Then, for Z ^ z\ the Schrodinger equation is more general than the Ince limit of DCHE. 
However, assuming a neutral target, we may form two pairs of solutions according to 



R^r) = -X^Ar) = e-^^/('^\^^-('/'^UU^ = r') (z = 1,2) 



(50) 



where on the right-hand side the Uiy represent the solutions with a phase parameter for the 
Ince limit of the DCHE, given in section 3.1. Then, taking into account that for this case 
q = k'^/A = fiE/{2h'^) and z = we find 



-i?i/{2r2)^- 



n=— oo 
oo 



[n + u + 



B, 



2n + 2p + 2: 



(51) 



V2 ^ 6«K2n+2.+i (±^fcr) ; 



^o^(^) ^ eB./(2.^V-2.-| g 5(2)('_:^y'vl/(n + z/ + 2-^,2n + 2z/ + 2;-^), 

n=-oo \ r J \ I r J 

(52) 

^oo(^) ^ giJ./(2.^)^-l/2 ^ (±^fcr) . 

?i=— oo 

From these expressions we obtain 

hm RZ{r) oc hm i?^(r) , -— < arg (izfcr) < — (53) 

where we have employed the limit (fT3|) for the modified Bessel functions. Thus, when r ^ oo, 
the solutions i?^ are bounded even if A; is a pure imaginary, since in this case exp(zA;r) or 
exp(— i/cr) goes to zero. At r = 0, Eq. (j39|l implies that 

lim<(r) - e-^^/^V^-(^/^ < arg^ < ^, 

^ r I ^^^^ 

limi?L(r) - e^^/^V(^/^)-^^ < arg (-^) < |. 
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Thence, if Bi is a positive real number, the first hmit goes to zero; if Bi is a negative real 
number, the second limit goes to zero. However, if Bi is a pure imaginary, we write 



Bi = iC, = 2 + ^ „ 



where C is real. Thus we find 

KWIoc|i?LWI~v^-o. 

Therefore, it is possible to find at least one pair of solutions for which both the solutions are 
bounded at the singularities. 



Potential without inverse sixth-power term. From Eq. (j^ we find that the substitutions 

z = r, x(r) = e-^^/(2rVi?2/2f^(^ ^ ^) ^-^^ B^ = -Afi{e'f , B2 = 2 

transform the Schrodinger equation ()49|) into 



dr"^ dr 



-1(1 + 1) - %iZ - z'Tee'r + + M4-6ao^;)(e^ 



i\2 



U = 0.(55) 



Thus, in absence of the inverse sixth-power term, the radial Schrodinger equation, even if we 
have a Coulomb term in the potential, may be solved by 



(56) 



where Uii,{z = r) are solutions with a phase parameter for the DCHE with z = r and B2 = 2 
(see Appendix B). As 



^ iZ - z'ree' k - ^ 



those solutions give 



n=— 00 
00 



kh' 



h 



R^^{r) = e^'^""-^ bn{T2ikr)''+'''^{n + u + l±ir],2n + 2u + 2;T2ikr)] 



(57) 



(5^ 



RUr) 



00 / R 



2r 



n=—oo 
00 



r 



n+v+l , 

^ I n + z/ + 1, 2n + 2z/ + 2; 

r 



J2 &n(T2ifcr)"+''^(n + z/ + 1 ± ir], 2ra + 2z/ + 2; T'^ikr). 



(59) 



Using Eq. (jSH)), we find 



lim R^^{r) oc lim i?^(r) ~ r^'"^ 



— ^ < arg (=Fzfcr) < — 



(60) 
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Thus, when r ^ oo, the solutions i?^ are bounded even if /c is a pure imaginary number, since 
in this case the behavior of exp{ikr) or exp(— zfcr) predominates over the other factor. At r = 0, 
by using Eq. (jH^ we get 

\imRUr)-e~^^/^^^-\ -^<arg^<^ 

2 r 2 ^g^^ 



, 2 -"°V r y ~ 2 
As Bi is a pure imaginary number, we find that 

KW|cxK(r)|~l 

Therefore, in this case we can form two pairs of solutions which are regular at the singular 
points, both pairs having the same series coefficients. For neutral targets {t] = 0) the previous 
results have already been found by Biihring who has treated the Schrodinger equation as a 
DCHE |16| l24j. Before this author, the Schrodinger equation (for neutral targets and an inverse 
fourth-power polarization potential) had been transformed into a Mathieu equation ^31 1^ . 
Thus, the Biihring approach is profitable since it works for ionized targets, too. In addition, as 
we have seen, for inverse sixth-power polarization potential, the Schrodinger equation may be 
transformed to the Ince limit of the DCHE, provided that the target is neutral. 



5. Final remarks 

We have constructed the differential equation (jH)) by applying the Ince limit, defined in Eq. (0), 
to a generalized spheroidal wave equation (GSWE). The Leaver limit (zq — >■ 0) of that equation 
has afforded Eq. (jSJ that turns out to be the Ince limit of a double- confluent Heun equation 
(DCHE) as well. The subnormal Thome behavior at 2 = 00, for the solutions of the these 
Ince limits of the GSWE and DCHE, distinguishes such equations from the original GSWE and 
DCHE hitherto considered in the literature. 

In section 2, a pair of solutions (with a phase parameter) for the Ince limit of the GSWE has 
been found as the Ince limit of a pair of solutions for the original GSWE. One solution is given 
by a series of hypergeometric functions and the other by a series of modifled Bessel functions 
of the second kind. Both solutions in that pair have the same series coefficients but different 
regions of convergence, as in solutions for the Mathieu equations. Other pair has followed from 
the first one by means of a transformation rule. Hence, four pairs of solutions without phase 
parameter have resulted from the truncation of the series with a phase parameter, that is, by 
restricting the summation index of the series to n > 0. 

In section 3, solutions for the Ince limit of the DCHE have been established by taking the 
Leaver limit of solutions for the Ince limit of the GSWE. These solutions, given by series of 
irregular confiuent hypergeometric functions and modified Bessel functions, present the appro- 
priate behavior at the irregular singularities z = and z = 00. Note, nonetheless, that in 
sections 3 and 4 we have dealt with expansions in series of modified Bessel functions only. 
Other possibilities may be investigated, specially solutions in series of Bessel function products, 
as these could have important properties as regards the convergence of the series. 

In the solutions without phase parameter for the Ince limits of the GSWE and DCHE, 
there are three possible forms to the recurrence relations for the series coefficients. This fact 
is relevant in itself and, in particular, is essential to recover solutions for the Mathieu equation 
from the ones for the Ince limit of the GSWE. 

The solutions we have obtained for the Mathieu equation are already known and exhibit 
the usual parity and periodicity properties. This includes also the solutions found by Poole, 
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given by two-sided series (—00 < n < 00) and having period 27rm, where m is any integer 
greater than 1. However, we note that other types of solutions for the Mathieu equations (and 
also for the Whittaker-Hill equations) are possible, since these equations may be considered as 
particular cases of both the GSWE and double-confluent Heun equations as well j2Z]. 

At last, notice that we have point out no application for Ince limit of the GSWE. Nev- 
ertheless, in section 4 we have seen that the Schrodinger equation ()49j) for the scattering of 
low-energy particles by polarizable targets leads to an DCHE and its Ince limit. The excep- 
tion is the Schrodinger equation with Coulomb and inverse sixth-power terms which requires 
solutions for a more general equation, possibly similar to an equation considered by Kurth and 
Schmidt in [28j. 

I thank Herman J. Mosquera Cuesta for his careful reading of this manuscript and valuable 
suggestions. I also thank the participants of the ICRA-BR "Pequenos Seminarios" for discus- 
sions and insight on potential physical applications of the results of the present investigation. 



Appendix A: Degenerate DCHEs 



Let us show that DCHE 



+ (5i + ^2^)^ + (^3 - 2r/cu^ + J'z') f/ = 0, (Si ^ 0, ^ 0), 



for 5i = and/or = degenerates into a confluent hypergeometric equation or an equation 
with constant coefficients. Thus, if _Bi = and 7^ 0, the substitutions 

y = -2iuz, U{z)=e-y/S''f{y). a" - {I - B^)a + = 
give the confluent hypergeometric equation 

,_^ + |(2a. + B.)-,l-^-(., + a+^)/ = 0. 

If i?i 7^ and u = 0, the change of variables 

y = Bjz, U{z)=y^g{y), /?2 _ (5^ - 1)/? + i^s = 

leads to 

d^q , ^do 

,_^+|(2;3 + 2-B,)-,|^-,a9 = 0. 

If Bi = uj = 0, we find an equation with constant coefficients by taking z = expy. 
Now let us show that the Ince limit of the DCHE 

(pTJ dU 

+ + ^2^):r + + qz)U = 0, (g ^ 0, 5i ^ 0) 
az"^ az 

also gives degenerate cases if g 7^ and/or Bi ^ 0. In fact, if g = and Bi 7^ 0, this equation 
is equivalent to the DCHE with u = and i?i 7^ 0. If g 7^ and Bi = 0, the substitutions 

e = ±2zv^, f/(z) = e'-^^T(0 

reduces the equation to the modified Bessel equation 



(1 - ^2)^ - 453 + e 



d^T dT 



2 ^ D I e2 



T = 0. 



de d^ 

Finally, for g = i?i = 0, we find again an equation with constant coefficients by taking z = exp y. 
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Appendix B: The solutions in series of Bessel functions 

The solution U^{z) in series of Bessel functions can also be constructed as follows. We perform 
the substitutions 

^ = ±2t^, Uiz)=e-^'YiO (Bl) 
in the Ince limit of the GSWE (jH)). This yields 

d'^Y 4g(2o - 25i - 2B2Z0) dY 



,d^Y dY 

h £ 

de di 



de 



+ 



4g(l - B,)^-^l±^^ + (1 - + Aqzo - AB, 



ds, 

Y. 



Now we expand Y[^) according to 

00 

^(0= E &i'^^A(0, A = 2n + 2z/ + l, 



(B2) 



(B3) 



where K\{^) denotes the modified Bessel function of the second kind 17^. The last equation 
and dHU) afford the solution U^{z). 

When we insert ()B3|) into ()B2|) . we use some difference-differential relations derived from 
the properties of Kx [T7j. Thus, we have 



on the left-hand side and 



de i d^ 

on the right-hand side. This gives 

2[l-2fi2-(2fii/^o)]" 



A- 1 



9^0 E 

ji=— 00 
00 

+ E 



1 + 



A- 1 



A^ + 453 -2g^o- (1-52)' 



2gzo[l-252-(2fii/^o)] 



n=— 00 

00 



A- 1 



n=— 00 

00 



QZo E 



1-2^2 



2Bi 



A + (1 - 52) 1 + ^2 



2Bi 



To remove the term 4/^^(0/^^ ^^e right-hand side we use the relation 



e A(A-l) (A-1)(A + 1) A(A + 1)' 
Then, reminding that A = 2n + 2z/ + 1, we find 



Y: ««i&i')i^2n+2.-l(0+ E /5i^^&i^)i^2„+2.+l(0+ E 7Sl&i'^^2n+2.+3(0 = 0, (B4) 
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where the coefficients and 7^^) are just the ones given in equations ()16b|) . To get the 

recurrence relations with the form given in piaj) . we change n — > m + 1 and n m — 1 in 
the first and third terms, respectively. After this, we equate to zero the coefficients of each 
independent K2m+2u+i{0- 

On the other hand, to study the convergence of the series, we apply a Perron-Kreuser 
theorem for the minimal solutions of the recurrence relations for fe^^-* and obtain (if zq 7^ 0) 



Using also the relation J7] 



lim £!l±i = lim ^ = (B5) 

Un Un 



lnni..(0 4r(A) f 



and ir_A(0 = i^A(0> we get 



Hence, we have 



i^'2n+2i.+l(0 ^2n+2;/+l(C) 



^n+l-^2n+2t/+3(0 _ &n-l-^2n+2!^-l (0 _ ^0 



n^oo 



2; 



Therefore, by the ratio test the series converges for \z\ > \zq\. In ()B5|) we have supposed that 
2^0 7^ but, if Zq = 0, we find 



hm — TTT- = hm — ttt- = 

j.^ ^i+l-^2n+2i/+3(0 _ j.^ ^i-l-^2n+2iy^l(0 _ 



n^oo 



On -^^2n+2;/+l (0 &5.'^i^2n+2.+l(0 



Thus, in this limit the series converges for |z| > and per se this result is already included in 

\z\ > \Zo\. 



Appendix C: Solutions for the DCHE of section 4 

The Leaver-type solutions for the DCHE Q present some simplifications for B2 = 2. The solu- 
tions given in Ref. |7] are expansions in series of regular and irregular confluent hypergeometric 
functions. However, to obtain the expected behavior at the singular points z = and z = 00, 
we have to choose the irregular functions. Then, by using the same notation of sections 2.1 and 
3.1, we find that for B2 = 2 the ffist pair of solutions with a phase parameter is given by 

^{'.(^) = e'"" E — ^U + z/ + l,2n + 2z/ + 2;— , 

n=-oo \ Z J \ Z J 

CI 

00 

U^{z) = e'^^ hn{-2iujzY+''m{n + u + l + iT],2n + 2u + 2]-2iujz), 
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and the second pair takes the form 

oo / j^^ V n+u+l , j^^ 



U^^{z) = e'^'+^ bni -) ^ (n + u + l,2n + 2u + 2; 

n=-oo \ Z J \ Z , 

U^{z) = e''^'+^ K{-2iujzf+'''^{n + u + 1 + ir],2n + 2u + 2; -2iujz). 

n=—oo 

Then, we see that the two pairs have the same series coefficients 6„ and the coefficients in the 
recurrence relations pia|l are simply 

«n = '^B, , /5„ = i?3 + (n + v){n + z/ + 1),7„ = iujB, (g^) . (C3) 

In these solutions v cannot be integer or half- integer and the [7°^ converge for any finite z, 
whereas the converge for \z\ > 0. Note, moreover, that the irregular confluent hypergeo- 
metric functions that appear in f/°, could be rewritten in terms of modified Bessel of the second 
kind by using the definition In the solutions the confluent hypergeometric functions 

could be rewritten in terms of the Hankel functions H^^^ but only if r/ = (neutral target, in 
the problem of section 4). For this we have to use the relation [1^ 



The asymptotic behaviors of the solutions given in (Cl-2) may be found by using Eq. (j39j) . 
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